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Method and Apparatus for Measuring the Profile of Small Repeating Lines 



This invention generally relates to a method of characterizing the size and shape of 
small periodic features on a substrate. More particularly, it relates to analyzing light 
reflected from a grating, or transmitted through the grating, to obtain line profile 
information. 



In microelectronics, accurate measurement of feature profiles (i.e. line width, line 
height, space between lines, and sidewall shape) are very important to optimizing device 
performance and chip yield. Measurements are needed at many steps in manufacturing to 
assure that critical dimensions, line profiles, and feature depths are under control. 
Historically, measurements have been accomplished with the following technologies: 

Optical imaging using the resolving power of optical microscopes and image 
processing to measure small features. However, features smaller than the resolving power 
of the microscope can not be measured, nor can the line profile. 

Electron-beam imaging, particularly the scanning electron microscope (SEM), 
greatly improves resolution over the optical microscope. However, like optical imaging, 
top-down SEM imaging does not provide profile information. While imaging a 
cross-sectioned samples does provide profile information, cross sectioning is destructive, 
costly, and labor intensive. In addition the signal processing used in SEM imaging 
introduces uncertainties. Furthermore, electron beams can damage the sample, and this is 
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especially the case when sensitive materials such as photoresist are imaged. In addition, 
electron beam charging can seriously distort measurement signals, causing beam position 
deviations. To overcome this problem, conductive coatings are commonly used but 
coatings or their removal can damage the sample or effect measurement accuracy. 

Scanning force microscopy (SFM) utilizes a very small mechanical probe to sense 
minute atomic forces that act very close to the surface of materials, and SFM can provide 
high resolution topographical information, including line profiles. Unfortunately, SFM is 
slow and operator-intensive, and the quantitative science of the probe tip is uncertain and 
unreliable. 

Surface-contact mechanical-probe technologies, such as surface profilometers, 
have been used to profile large structures but do not offer the resolution or sensitivity 
required by today's characterization requirements. Furthermore, contact-probes distort 
surfaces. In addition, mechanical stability requirements prohibit the use of probes small 
enough to accommodate the submicrometer sizes commonly measured. 

To measure the width of a line the location of each edge of its profile must be 
defined, and an arbitrary, qualitative edge detection model is usually used in each of the 
above measuring techniques. While this arbitrary measurement point is typically calibrated 
to a cross-section, manufacturing process changes and normal manufacturing process drift 
can invalidate the edge model and introduce significant measurement error. 

Furthermore, none of the above techniques can characterize planarized, buried 
sub-micron structures that commonly are found in microelectronic or micro-machining 
manufacturing processes. Features may be formed that are later covered with a transparent 
(or semi-transparent) layer of material which is then planarized. If the feature sizes are 
below the resolution limits of optical techniques and the flat top surface prevents SEM and 
AFM imaging, none of the above well known techniques will work. 
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Two optical techniques are available for planarized buried structures. 
Scatterometry has been used for characterizing periodic topographic structures in which a 
light beam (typically a laser) illuminates an area to be characterized and the angular 
distribution of elastically scattered light is used to measure line widths on photomask and 
silicon wafer gratings. However, scatterometry requires a high degree of mechanical 
stability and accurate measurement of the laser beam impingement angle on the grating, 
and these mechanical limitations contribute to the uncertainty of scatterometry results. 

Confocal laser imaging is another optical technique that permits characterization of 
planarized buried structures. Focus on the sample is progressively adjusted with a 
monochromatic, high numerical aperture (NA) optical imaging system, and those focus 
adjustments provide the line profile. The high NA lens has a very small depth of focus, and 
therefore, a physical profile can be obtained by moving through a series of focus settings 
up or down the features, and recording those height variations. However, the technique is 
limited by the resolution of the laser beam, or the size of the laser spot on the sample. 
Furthermore, the use of a high NA objective lens, providing a large angle to the incident 
beam limits the ability to characterize high aspect ratio features that are commonly found 
in semiconductor manufacturing. 

Thus, a better solution for characterizing microstructure geometries is needed to 
determine the two dimensional physical profile of a line with a method that is rapid, 
accurate, reproducible, and reliable, and this solution is provided by the present invention. 



It is therefore an object of the present invention to provide an optical method of 
determining a line profile of a sub-micron structure. 



Summary of the Invention 
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It is a feature of the present invention that the line profile is determined by 
comparing measured diffraction radiation intensity from a periodic structure on a substrate 
with modeled radiation intensity predicted from a theoretical line profile. 

It is a feature of the present invention that line profile parameters of the model are 
varied until the predicted intensity converges to the measured intensity. 

It is a feature of the present invention that a rigorous and detailed two dimensional 
analytical solution is obtained for the line profile. 

It is an advantage of the present invention that the line profile is reliably, 
reproducibly, rapidly, and nondestructively, obtained. 

These and other objects, features, and advantages of the invention are 
accomplished by a method of determining the profile of a line comprising the steps of: 



(a) providing a substrate having a repeating structure comprising a 
plurality of lines, said lines having substantially identical profiles; 

(b) illuminating said repeating structure with radiation wherein said 
radiation diffracts, said diffracted radiation having an intensity; 

(c) measuring said intensity; 

(d) providing a model structure on a data processing machine, said 
model structure comprising a repeating structure on said substrate, 
said model structure comprising a model profile; 

(e) mathematically predicting a predicted diffracted radiation intensity 
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when said model structure is illuminated with said radiation; and 
(f) comparing said predicted intensity with said measured intensity. 

Brief Description of the Figures 

The foregoing and other objects, features, and advantages of the invention will be 
5 apparent from the following detailed description of the invention and accompanying 
drawings, in which: 

FIG. Lis a flow chart of the steps of the present invention, showing the iterative 
process for calculating a line profile; 

FIG. 2-is an ideal grating having a vertical line profile; 

FIG. 3 is a an intensity versus wavelength curve for the ideal grating of FIG. 2; 

FIG. 4 is a an ideal grating having a sloped line profile; 

FIG. 5 is an intensity versus wavelength curve for the ideal grating of FIG. 4; 

FIG. 6 is an overlay of the intensity versus wavelength curves shown in FIGS. 3 
and 5 illustrating the effect of a change in grating profile on the intensity versus 
15 wavelength curve; 

FIG. 7 is a cross section illustrating how a line profile is approximated by a stack 
of slabs, each slab having a width, a height, and an index of refraction, and showing the 
grating periodicity; 
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FIG. 8a is a seed line profile having an S shape corresponding to the data of Tables 
1 and 2; 

FIGS 8b and 8c are top and bottom portions of the line profile of FIG. 8a; 

FIG. 8d is the seed profile of FIG 8a after scaling the data of Tables 1 and 2. 

FIG. 9 is a schematic of the apparatus of the present invention; 

FIG. 10 illustrates diffracted angle from a grating as a function of incident angle 
and grating periodicity; 

~ , FIG. ¥$£APL code for calculating the line profile of the present invention; and 

w Xa<o * 

fp 

^ FIG. 12 is a plot of measured reflectivity and theoretical reflectivity (TM and TE) 

^0 10 and the corresponding line profiles for TM and TE calculated according to the method of 

ru 

gO the present invention. 
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FIG. 13a is a top view of a two dimensional grating of trenches, such as DRAM 
capacitor trenches for which two dimensional line profiles can be calculated according to 
the present invention. 

15 FIGS. 13b, 13c are orthogogal cross sections of the grating of FIG. 13a. 

FIG. 14a is a cross sectional view of a semiconductor substrate having a doping 
that varies with depth into the substrate, the depth profile of which can be calculated 
according to the present invention. 

FIG. 14b is a cross sectional view of a substrate having film layers of varying 
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thickness and index of refraction, the depth profile of which can be calculated according to 
the present invention. 

FIG. 14c shows the real and imaginary parts of the index of refraction versus depth 
profile calculated according to the present invention for the doping profile of FIG. 14a. 

FIG. 14d shows the concentration versus depth profile as converted from the index 
of refraction curves of FIG. 14c. 



The present invention is a method for nondestructively determining the 
topographical cross-section of lines on a substrate which provides line thickness, line 
width, and the shape of the line edge (the line profile). While a repeating structure, or 
grating, is required for the measurement, the method uses broad band illumination, does 
not involve contact with the substrate and can equally be used for buried planarized 
gratings. The method takes advantage of available parallel processing computer 
capabilities for providing rapid line profiles. 

Diffraction gratings are often displayed with perfectly square, sine wave, sawtooth, 
triangular, or other ideal edges. But, in fact, the edges of all gratings deviate from the 
ideal, often significantly. There are important applications in semiconductor manufacturing 
where knowledge of the shape of the line edge, the line profile, is useful for monitoring or 
controlling a process. It is therefore one primary goal of this invention to provide an 
improved method of determining this line profile. In addition, there are other applications 
in semiconductor manufacturing where composition profile versus depth, such as doping 
concentration versus depth is needed, and the invention also satisfies this need. 




Detailed Description of the Invention 
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The flow chart of FIG. 1 illustrates the method of the present invention. In the first 
step (step 100), convergence criteria are established between measured and calculated 
curves acquired in steps described hereinbelow. Next a substrate having a repeating 
structure, such as a grating, is illuminated with broad band radiation. Diffracted radiation 
5 is collected, measured, normalized to the incident radiation, and recorded as a function of 
wavelength to provide an intensity versus wavelength curve (step 102). Then, an initial 
model (or seed model) of the line profile of the grating, a model of the broad band 
radiation shined on the grating, and a model of the interaction of the radiation shined with 
the model grating is provided to a data processing machine (step 104). The data 
1 0 processing machine uses Maxwell's equations to calculate a model diffracted intensity 
versus wavelength curve (step 106), and the measured intensity curve is then compared 
with this modeled intensity versus wavelength curve (step 108). If agreement between the 

□ curves within the convergence criteria of step 100 is not found, the line profile in the 

C8 

s « model is then adjusted (step 1 10) and the model intensity curve recalculated to attempt to 

f 3 15 improve agreement between the measured and calculated intensity curves (steps 1 06 and 

l D 108 repeated). The model is repeatedly adjusted and the intensity recalculated until 

yjj agreement, within the convergence criteria established in step 100, between the two 

f intensity versus wavelength curves is achieved. The theoretical profile is the actual profile 

□ to an accuracy determined by the extent of the prescribed convergence limits, 

W 

y\ 20 measurement accuracy, and the extent to which the seed, as modified by scale factors, can 
^ approximate the actual profile (step 112). 

The invention takes advantage of the fact that the magnitude of the energy 
diffracted into the various diffraction orders is strongly influenced by the line profile in 
addition to line width, line spacing, and thickness and optical properties of the materials 
25 comprising the grating. Because of the sensitivity of the diffraction pattern to these 

additional parameters, the present invention provides a means to extract one or more of 
these parameters in addition to the line profile. 
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The invention extracts a profile-dependent signal from a grating by measuring the 
intensity of diffracted radiation over a small solid angle. Thus, only a portion of the 
diffracted radiation is measured, such as the radiation in the zeroth order of the diffraction 
pattern. Because details of the line profile affects the distribution of energy among the 
5 diffraction orders, the line profile may be obtained by measuring the intensity in only one 
of the orders. 
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FIG. 2 shows a diffraction grating with idealized vertical edges. The interaction of 
this grating with broad-band illumination (light having a range of colors) can be modeled 
using Maxwell's equations, and the reflected intensity as a function of wavelength for the 
10 zeroth order (the order near normal to the surface) can be calculated, as shown in FIG. 3. 
Preferably the grating has at least 5 lines. 

FIGS. 4-6 illustrate how changes in the line profile can change the energy so 
diffracted. In FIG. 4, a grating similar to that of FIG. 1 is used, but with line edges having 
a slope of about 85 degrees. The intensity versus wavelength curve predicted by 



ry 

ijQ 15 Maxwell' s equations for the grating of FIG. 4 is shown in Fig 5 . To further illustrate the 



difference in intensity versus wavelength curves due to the change in line profile, the 
curves of FIGS. 3 and 5 are superimposed in FIG. 6. The marked difference in intensity 
versus wavelength demonstrated in FIG. 6 is taken advantage of in this invention to 
determine line edge profiles. 

20 In the model of the line profile used for predicting intensity versus wavelength, line 

profile information is expressed as a set of stacked slabs of material, as shown in FIG. 7. 
Each slab is defined by a width w, , a height d, , and an index of refraction n ( . Any profile 
can be approximated by a series of such stacked slabs. The profile can be made arbitrarily 
smooth by including a sufficient number of slabs. However, calculation efficiency 

25 demands that a minimum number of slabs be used, and by properly adjusting the height 

and width of each slab, any arbitrary profile can be represented with a minimum number of 
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Thus, the method begins with an initial model of the line profile expressed 
as a series of slabs. This profile is used as an input to a computer program that predicts the 
percentage of reflected energy diffracted into the zeroth order over a range of wavelengths 
for the model profile. The predicted intensity is compared to the measured intensity as a 
function of wavelength (normalized for the incident intensity) and the slab widths and 
heights are adjusted until agreement between the predicted and measured intensity versus 
wavelength curves is achieved. The final result of the analysis is a stacked set of slabs that 
represents the line profile of each line of the grating. 

The inventors of the present invention have found that only about 20 slabs are 
needed to adequately represent many profiles found on semiconductor wafers. Since each 
slab has two independent values associated with it (a width and a height), two numbers are 
needed to describe each slab, and about 40 independent variables are needed to 
approximate a profile. While, in principle, these 40 variables could be adjusted as 
described above to yield best agreement between measured and predicted diffracted 
intensity, the computational task of determining the edge profile has been found to be 
smaller by reducing the number of variables. 

In the preferred embodiment of the calculation, the inventors have reduced the 
number of variables by (1) dividing the model line profile into two or more sub-profiles; 
and (2) providing a numerical model of each sub-profile wherein a relatively few scaling 
factors are used to adjust all slab widths and heights within a single sub-profile 
simultaneously. Thus, the problem is reduced from 40 variables to one having a few 
variables. 

For example, a single scaling factor is used for adjusting all of the width points 
within each sub-profile and another scaling factor is used for adjusting all of the height 



BU9-96-041 



10 



points of the "S" line profile of FIGS. 8a-8c. This "S" profile (FIG. 8a) can be viewed as 
having an upper curved region (FIG. 8b) and a lower curved region (FIG. 8c), each with 
different curvature. The curvature of the bottom of the S and the top of the S are not 
known in advance and they may be quite different in magnitudes. This S profile is divided 
into 2 sub-profiles each formed of many slabs. The widths of all slabs in each sub-profile 
are varied according to a scaling factor that permits a wide range of sub-profile shapes. 
The two scaling factors for the two sub-profiles are allowed to vary independently. The 
height of all slabs in each sub-profile is adjusted by a third scaling factor, and this same 
scaling factor can be used for both sub-profiles. Thus, the entire S shaped edge profile can 
be described by three independent variables. 

Table 1 illustrates an initial estimate (seed) of the basic shape of the upper 
sub-profile of the "S "-shaped profile as a table of layer index, X-position, and slab 
thicknesses (d,). The X-position defines the edge shape and has a reference to zero for 
both top and bottom such that sub-profile shapes may be modified by scaling factors as 
described below while maintaining continuity between the sub-profiles. 



X d 



1 


0.8 


0.00025 


2 


0.75 


0.00025 


3 


0.7 


0.00075 


4 


0.65 


0.00125 


5 


0.6 


0.002 


6 


0.55 


0.00325 


7 


0.5 


0.00475 


8 


0.45 


0.007 


9 


0.4 


0.0095 


10 


0.35 


0.013 
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11 


0.3 


0.01725 


12 


0.25 


0.02275 


13 


0.2 


0.029 


14 


0.15 


0.0365 


15 


0.1 


0.046 


16 


0.05 


0.0565 


17 


0 


0.5 



Table 2 is a sub-profile seed for the lower portion of the "S"-shape. 
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26 


-0.45 


0.007 






27 


-0.5 


0.00475 






28 


-0.55 


0.00325 






29 


-0.6 


0.002 






30 


-0.65 


0.00125 




25 


31 


-0.7 


0.00075 






L=32 


-0.75 


0.00025 



The upper and lower seed profiles of Tables 1 and 2 are shown in FIG. 8a. The 
BU9-96-041 12 



magnitude of curvature of the upper sub-profile (Table 1) is changed by allowing the X 
axis to change by a linear scale factor, so that the X variation of a 'new' sub-profile is 
given by 



where X, is the initial value of X for the fth slab, m, is the scale factor, and X ( ' is the value 
of X obtained for the 1th slab after applying the scale factor. Slab thickness values d, in 
Tables 1 and 2 are retained while the X values change according to the m : scale factor. 



In a similar manner, the lower sub-profile can be adjusted by allowing the X values 
associated with it to change by a different linear scale factor 



where X ( is again the initial value of X for the fth slab, m 2 is the scale factor for the lower 
sub-profile, X,' is the value of X obtained for the fth slab after applying the second scale 
factor, m 2 , and L is the index number of the bottom slab, in this case slab 32. 

One additional scale factor is recommended for the slab thickness values such that, 



The scale factors in the method given above allow for a wide range of 'S' shaped 
curves to represent the line profile that are formed from a seed of initial X and thickness 
values as modified by three scaling variables, m u m 2 , m 3 . FIG. 8d shows the modeled line 
profile for m t = 2, m 2 = 0.01, and m 3 = 1.2, illustrating the effect of the three scale 
factors on the seed curve. In addition to linear scale factors, higher order scale factors or 
scale factors having other functional forms can be used to modify the seed profile. 



X/ - m, X, 



f = 1 ... 17 



X,' — m 2 Xj ( — 1 8 ... L 



d|' = m 3 d, I = 1 . . . L 
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In practice an additional parameter W will be required to set the effective line 
width of our grating elements and to define the separation between the symmetric left and 
right edge profiles of a single line of the grating. The parameter W will explicitly define the 
width of each slab when used in conjunction with the scaled edge profiles. For example, 
the width of each slab for the lower profile will be given by 



w ( = 2X, + W l=\...L 
With this additional parameter, we have reduced the number of parameters to be varied 
from about 40 to 4. The optimization program is now used to successively change these 4 
profile parameters until the gap between measured intensity and predicted intensity versus 
wavelength is minimized. The inventors have found that the seed chosen to initiate the 
iterative process is important for the success of the iterative procedure. Thus, an 
approximate shape obtained from a cross section or atomic force microscopy profile that 
approximates the shape of the profile will improve the chance of success and decrease the 
number of iterations required. 



Fig. 9 shows the apparatus used to measure the actual diffraction vs. wavelength 
curve of a particular grating. Broadband illumination source 10 is projected on line 
diffraction grating 12 on substrate 13 through optical apparatus 14, such as partially 
silvered mirror 16 and optical microscope objective 18. Light incident on grating 12 from 
objective J^is diffracted into many diffraction orders, of which orders 0, 1, and 2 are 
shown for reflected and transmitted diffracted light. The analysis of diffracted energy 
reflected from absorbing substrates, such as silicon substrates, will be described in detail 
here. In the present invention, the diffracted light energy associated with one or more of 
the low orders is collected by microscope objective 18, transmitted through partially 
silvered mirror 16,-potarizedi20, and collected by spectral separator 22a and detector 22b. 
Polarizer 20 provides a single polarization •fteajythe collected light energy to detector 22b, 
that polarization preferably being either transverse-electric (TE) or transverse-magnetic 
(TM). While polarizer 20 is shown in detector light path 24, polarizer 20 can also be 
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placed in source light path 26. Although polarized light is preferred, unpolarized light or 
other polarization states are acceptable if the analysis can account for the situation. 

Within spectral separator 22a, the polarized light energy is separated into its 
spectral components by a monochromator (not shown) and the intensity as a function of 
wavelength is measured by detector 22b, both of which are well known in the art. The 
relative reflectivity at each wavelength is then computed from the detected intensity as a 
function of wavelength. The computation involves a correction for the intensity of the 
incident light at each wavelength. Since the illumination source does not produce a " 
constant intensity at each wavelength and since the optics of the system change the 
intensity of different wavelengths differently, it is necessary to normalize the system. This 
is achieved with a reflectivity reference, such as polished silicon. The reflected intensity as 
a function of wavelength from the reflectivity reference mounted in the system is measured 
and the signal compared with the known reflectivity of the reflectivity reference as a 
function of wavelength. A normalization factor is thereby generated at each wavelength 
which can be applied to correct the intensity of signals from an actual grating where the 
reflectivity is to be measured. 

Incoming light to grating 12 and diffracted light reflected from grating 12 are 
shown in FIG. 10. The incident light having wavelength A ; arrives at an incident angle 0; 
from normal. The diffracted light is reflected at angle 0 S which is given by 

6 s = sin "'(sine,. - 

where m is an integer referred to as the "diffraction order" and A is the grating period. The 
diffracted light reflects with a wavelength unchanged by the interaction. 

The present inventors have found that choosing an optical path wherein the 
incident light is normal to the diffraction grating simplifies the calculation of line profile. 
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They have found that the calculation is also tolerant to a small range of angles symmetric 
around such a normal incident beam. Such illumination is conveniently achieved by using a 
low numerical aperture objective in the optical microscope. In practice, a numerical 
aperture less than .08 is used. 

In addition, the present inventors have designed optical apparatus 14 such that a 
significant percentage of the diffracted energy is diffracted at angles too large to be 
collected by objective lens 18; that is, the diffraction angle of the reflected light is outside 
the numerical aperture of the lens. From the diffraction equation above, it can be seen that 
the zeroth diffracted order is scattered at the same angle as the incident light. In the 
preferred embodiment, objective lens 18 and the wavelength range are selected so that 
only the zeroth diffraction order is collected by the apparatus. 

A line profile modeled as a stack of slabs, each with a specific width, height, and 
complex refractive index used in the present invention, illustrated in FIG. 7, provides a line 
profile resolution determined by the number of slabs. By increasing the number of slabs 
comprising the line the shape of the edge profile can be made arbitrarily smooth, while 
decreasing the number of slabs decreases the time for computing the line profile. 

TM polarization solution 

In the APL code contained in FIG. 11, Maxwell's equations are solved to provide 
the electric and magnetic fields and light intensity reflected by a grating in each order for 
the optical system illustrated in FIG. 10 using a technique generally described in a paper 
by Morham in the J.Opt.Soc.Am., Vol.12., No.5, May 1995. The reflected diffraction 
efficiencies (the percentage of light reflected for each diffraction order or reflectivity) are 
also calculated in the APL code. This code provides the partial solution of Maxwell's 
equations for reflected diffraction, as needed for line profiles on substrates such as non- 
transparent semiconductor wafers. While many methods for solving Maxwell's equations 
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are known, the full solution, including both reflected and transmitted light, can be derived 
following the description in that paper and others by the same author. 

In the code of FIG. 11 permitivity harmonics (E) are calculated from the Fourier 
harmonics of the input parameters (n, f) as shown in the function PERMITIVITY and 
5 PERMPRIME. Wavenumber equations (K) are derived in matrix form in the function 
WAVENUMBER from the input wavelength, grating period, and angle of incidence. 
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The coupled wave equations, in matrix form, (B; COUPLEDWAVE line 28) result 
from the requirement that electromagnetic fields satisfy Maxwell's equations in the grating 
layer. EIGENSTUFF performs the eigenspace calculations to solve the coupled wave 
10 equations, producing the eigenvector diagonal matrix (W) and the square root (positive 

real part) of the eigenvalues as the diagonal matrix (Q) combining with the permittivity (E) 
to produce the product matrix (V; COUPLEDWAVE line 32). 



S A system of 2n(L+l) equations (where n is the number of space harmonics 

yg retained in the solution and L is the number of grating layers used to form the profile) 

j\, 15 results from the matching of electromagnetic fields at the grating layer boundaries. 



The standard transmittance matrix solution approach can become unstable due to 
finite computing precision of matrix inversion. Therefore, the solution employed in the 
code uses the numerically stable transmittance approach to determine the diffracted 
reflected amplitudes only. This matrix solution method utilizes a set of matrix calculations 
20 that cascaded from the first layer to the layer (function FANDG called from 

COUPLEDWAVE line 39). The resulting wave matrix (R; COUPLEDWAVE line 40) is 
used to calculate the diffracted efficiencies (DERTM; line 42). Calculations for the 
transverse electric field follow similarly in lines 47 through 60. 

The function COUPLEDWAVE is invoked by issuing the APL command, 
BU9-96-041 17 
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COUPLEDWAVE WL, where WL is a required argument; its value being the wavelength 
at which to evaluate the theoretical profile. COUPLEDWAVE, as configured in FIG. 11, 
is set up to compute TM diffraction. To change to TE, remove the comment symbol from 
line 26 & 63 and add a comment symbol to line 62. COUPLEDWAVE also requires 
several other variables to be defined in the work space and these are described in Table 3: 



Table 3 



Variable Name 


Format 


Purpose 


LAYER 


An N by 3 matrix; where N 


Describes the theoretical 




is the number of layers used 


profile of the grating 




to duiiq tne tneoreticai 






profile, column 1 is the 






layer film ID (see 






CAUCHY), column 2 is the 






layer width (nm), and 






column 3 is the layer 






thickness (nm). 




1 IIS VJJT JGIVl \JU 


r\. sCdlal II1II1J. 


1UC J^IaUU)4 pcllUU 


ORDERS 


A scalar integer. 


The number of positive 
diffracted orders retained in 
the analysis. 


TH 


A scalar (degrees). 


The angle of incidence. 


WL 


A scalar (nm). 


The incident wavelength. 


CAUCHY 


An M by 4 matrix where M 


A lookup matrix of Cauchy 




is any number of film types 


coefs. to calculate film 




for which Cauchy coef. 


refractive index (user 




data is known; column 1 is 


supplied). 




the film ID, cols. 2-4 are 






Cauchy coefs. CI, C2, & 






C3. 





BU9-96-041 



18 





SI 



Is a P by 3 matrix where P 



The complex refractive 



is the number of refractive 



index of silicon used as the 



index values included, 
column 1 is the wavelength, 



substrate (user supplied). 



cols. 2 and 3 are the real 



and imaginary parts of the 
complex refractive index at 
the given wavelengths. 



GRATING PROFILE 

FIG. 12 shows a line profile result of the present invention for reflectivity versus 
wavelength data taken from a photoresist line on a silicon substrate. The profile was 
calculated independently from two normal polarizations, transverse magnetic (TM) and 
transverse electric (TE). In the case of TM, the magnetic field vectors are perpendicular to 
the lines of the grating whereas electric filed vectors are always perpendicular to magnetic 
field vectors. Therefore, TE polarization has electric field vectors perpendicular to the 
lines of the grating. The different polarization orientations may be obtained by rotating the 
polarizer in the optical path so only light having the polarization of interest is collected. 

To minimize analysis time, it is convenient to select either the TE or the TM 
polarization. In this way only the normalized reflected intensity for a single polarization 
needs to be calculated. In FIG. 12 the line profile as calculated with both TE and TM is 
shown. The figure shows close agreement between the calculated intensity versus 
wavelength curves and the actual intensity versus wavelength curves for both TE and TM. 
Both give almost identical line profiles. 

The present invention can also be applied to calculating the line profile of trenches. 
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For example, the edge shape profile of a phase shift structure on a quartz mask, which is a 
trench etched in the quartz, can be calculated with the present invention. In this case, 
either the transmitted or the reflected diffracted energy can be used to determine the depth 
and edge profile. Similarly, edge profiles of trenches in a semiconductor substrate can be 
5 calculated. Two dimensional line profiles of gratings, such as contacts or DRAM cell 
capacitor trenches, as shown in FIG. 13a- 13c can also be calculated by extending the 
formalism for solving Maxwell's equations to two dimensions. 



In addition to calculating line profiles, the present invention can be applied to 
calculating other types of optically sensitive profiles, including doping profiles and 
10 composition profiles. In semiconductor manufacturing, it is important to control the 

concentration of dopants as a function of depth. FIG. 14a an implanted semiconductor 

O having a doping concentration c(z) near the surface that varies with depth, z. It is well 

B 

j3 known that the optical characteristics of silicon vary significantly with dopant 



^ concentration in the infra-red optical range, from a wavelength of about lum to a 

€! 1 5 wavelength of about 40 urn. In the first step a grating of suitable period is etched into the 

ru 

yp doped silicon to diffract the first order outside the range of the detected angles, as 

j 3 . described hereinabove. Preferably the sidewalls of the grating have a straight vertical line 

O profile, as shown in FIG. 14b, to simplify the calculation exclusively to the varying optical 

\1 properties. In the next step, the grating is illuminated with the broad band infrared 

rj 20 radiation. The diffracted normalized intensity versus wavelength is measured as described 
hereinabove. A seed model of the index of refraction depth profile is provided to a data 
processing machine and the process described hereinabove is used to calculate a 
normalized intensity versus wavelength curve for comparison with the measured curve. 
The model is adjusted to improve agreement and the calculation repeated until 
25 convergence is achieved as described hereinabove. The index of refraction depth profile of 
FIG. 14c is then converted into a doping profile using, for example, published tables 
relating index and doping concentration, as shown in FIG. 14d. 
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The present method can be applied to calculate the thickness and index of 
refraction of a stack of optical coatings, such as antireflective coatings on a lens. An edge 
of the lens or a monitor lens could be used for the grating. The present invention can also 
be applied to calculating the line profile and index profiles of buried structures, such as 
those shown in FIG. 15. The various indices and thicknesses of the layers can be 
determined as described hereinabove. 

While several embodiments of the invention, together with modifications thereof, 
have been described in detail herein and illustrated in the accompanying drawings, it will 
be evident that various further modifications are possible without departing from the scope 
of the invention. For example, a library of line profile curves can be used instead of or in 
combination with iterating from a seed profile. Nothing in the above specification is 
intended to limit the invention more narrowly than the appended claims. The examples 
given are intended only to be illustrative rather than exclusive. 

What is claimed is: 
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